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Abstract.  The  problem  of  a  two-phase  dispersed  medium  is  studied  within  the  kinetic  theory.  In  the 
case  of  small  and  undeformable  spheres  that  have  all  the  same  radius,  a  model  for  the  collision  operator 
of  the  Boltzmann  equation  is  proposed.  The  collisions  axe  supposed  instantaneous,  binary  and  inelastic. 
The  obtained  collision  operator  allows,  to  prove  the  existence  of  an  H  theorem  in  several  configurations 
according  to  the  assumptions  made  about  the  particles  and  particularly  in  the  case  of  a  diluted  suspension 
of  weakly  inelastic  collisions.  Because  of  the  complexity  of  the  non  linear  structure  of  the  collision  integral, 
the  Boltzmann  equation  is  very  difficult  to  solve  and  to  analyse.  It  is  therefore  interesting  to  introduce  a 
BGK  model  equation  to  study  qualitatively  its  solution.  In  order  to  be  physically  realistic  and  consistent 
with  the  assumptions  related  to  the  collisions,  a  collision  frequency  depending  on  the  particles  velocities  is 
chosen.  Moreover,  the  collision  frequency  is  expected  to  vary  strongly  with  the  particles  velocities.  Taking 
this  into  account,  a  BGK  model  is  written.  All  the  basic  properties  of  the  original  operator  are  retained. 


I  BOLTZMANN  EQUATION 

The  problem  of  a  two-phase  dispersed  medium  is  studied  within  the  kinetic  theory.  In  the  case  of  small  and 
undeformable  spheres  that  have  all  the  same  radius  {D/2),  a  model  for  the  collision  operator  of  the  Boltzmann 
equation  is  proposed.  Firstly,  we  consider  the  binary,  inelastic  and  instantaneous  collision  of  two  identical  hard 
spheres  of  radius  0/2  {Pi  and  P^).  Before  the  collision  we  suppose  that  the  particle  Px  is  centered  at  position 
xi  and  has  the  velocity  The  particle  P2  is  centered  at  position  and  has  the  velocity  ^2-  So,  the  relative 
velocity  is  p  =  ^2  -  and  the  impact  vector  is  fc  =  {xi  -  x^)/  ||  x\  -  £2  ||-  After  the  collision,  we  suppose  that 
the  particle  Pi  is  centered  at  position  xi  and  has  the  velocity  '.  The  particle  P2  is  centered  at  position  X2 
and  has  the  velocity  ^2  ^  and  the  relative  velocity  after  collision  is  g'  =  ^2'  —  ^1' ■  Furthermore,  because  of  the 
inelasticity  of  the  collision,  we  introduce  the  coefficient  of  restitution  e.  It  is  defined  by  the  following  relation 
between  the  components  of  the  relative  velocities  normal  to  the  plane  of  contact. 

{9  ■  k)  = -e{g' ■  k)  (1) 

If  e  =  1  the  energy  is  conserved  during  the  choc  and  if  e  <  1  energy  is  dissipated  during  a  collision.  By 
considering  the  momentum  balance  and  the  previous  assumption,  we  can  express  the  particle  velocities  just 
after  the  collision  in  term  of  those  just  before  and  reciprocally  ; 


(5  •  fc)  fc] 

(2) 

a 

(5 '  •  fc)  fc] 

(4) 

{g-  fc)  k 

(3) 

a 

[9'  ■  fc)  fc] 

(5) 

We  consider  the  following  Boltzmann  equation  (6)  for  a  suspension. 
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Ot  dxi 


where  F  is  the  external  force  (due  particularly  to  the  carrying  fluid)  acting  on  a  particle  and  where  we  have 
introduced  the  distribution  fonction  f{^,x,t)  which  is  defined  so  that  f{^,x,t)d^  is  the  probable  number  of 
particules  in  x,  at  time  t  with  velocities  in  the  volume  element  centered  in  ^  and  denoted  by  {^,d^).  We 
introduce  also  the  distribution  function  which  characterizes  the  statistic  of  binary  collisions  and  depends 
on  the  velocities  and  positions  of  two  particules  and  of  time.  It  is  defined  so  that 


,x{,^2,X2,t)d^ld£,2dxidX2 


(7) 


is  the  probable  number  of  pairs  of  particles  which,  at  time  t,  are  in  {x{,dxi)  and  {x2,dx2)  with  velocities 
respectively  in  (^i,d^i)  and  i^2,d^2)-  We  have  to  determine  C{f,f)  in  term  of  Using  a  very  classical 
method  (  [4])  in  kinetic  theory  based  on  the  change  in  {xi,dFi)  of  the  number  of  particles  with  velocity  in 
during  dt,  we  obtain  (  [5])  the  following  expression  for  the  collision  operator  of  the  Boltzmann  equation 


ot  ox  I 


(g-k)  +  Dk,t) 

-  Dk,t)'^  dkd(2 


e-"  Jh9  k)>0 


(8) 


We  choose  for  the  formulation  given  by  [4]  : 

+  Dk,t)  =  go{x{  +Dfc/2)/(^i,ari)/(^2,a;l  +  Dk)  (9) 

It  allows  to  take  into  account  the  dimension  of  the  particles.  This  expression  introduces  the  function  go{x) 
which  expresses  the  influence  of  the  density  of  the  cloud  of  particles  upon  the  probability  of  collision  :  go{x)  =  1 
for  a  diluted  suspension  and  go{x)  becomes  infinite  when  the  density  is  so  great  that  the  suspension  behaves 
like  a  solid.  Particular  cases  allow  the  validation  of  the  operator.  With  go{x)  =  1,  e  =  1  and  xi  =  X2  (punctual 
collisions)  the  collision  operator  of  the  usual  kinetic  theory  ( [3],  [4],  [10])  is  found.  If  we  suppose  that  goix)  =  1 
and  =  X2  the  operator  of  [14]  is  obtained.  At  last,  with  e  =  1,  the  Enskog  equation  (  [4]  and  [2])  is  obtained. 

The  function  i/'(^i,aii)  is  introduced  and  the  equation  of  moments  is  written  in  the  classical  way  of  the  kinetic 

theory.  We  put  ip{^i,Xi)  =  -ipi  with  i  =  1,2  and  ip{^[,x\)  =  with  i  =  1,2.  It  has  been  seen  in  [5]  that  the 
contribution  of  collisions  C(t/>)  to  the  transport  equation  can  be  written  in  the  same  form  as  the  expression 
directly  postulated  by  [9]  and  [4]  (for  dense  gas)  without  introducing  the  Boltzmann  equation  : 

^  [[  ig-k)\{rpi' -rpi)  f^^H^i,x"i,^2,x"i  -  Dk,t) 

+  -'<p2)  f^^\^2,xi,^i,xi  +£>fc,f)j  dkd^2d^i  (10) 

with  (2)  and  (3).  This  expression  provide  the  conservation  of  mass  (C(l)  =  0).  If  we  consider,  in  addition  to 
the  initial  assumptions,  that  the  diameter  of  the  particles  is  very  small  compared  to  the  characteristic  length 
(L)  of  the  flow,  and  thanks  to  an  asymptotic  expansion  in  terms  of  the  small  parameter  D  jL  «  1,  we  obtain 
a  relation  that  allows  to  prove  the  conservation  of  momentum  :  the  collision  term  of  the  momentum  balance 
equation  may  be  written  under  the  form  of  a  divergence  term  (  [9]).  Assuming  that  e  is  nearly  one  (weakly 
inelastic  collisions)  and  using  a  linear  combination,  we  can  express  C(^)  in  a  rather  symetrical  form  [5]  : 

^  2(1 +  e^)  [{^'i  ~  (^f^Hii',xi,i2',xi+ Dk,t) 

-  f^^H{i,xi,{2,xi  -Dk,t)^  +(V'2'-V'2)  -  Dk,t) 

-  f^^\^2,xi,^i,xi  +  Dk,t)"j^dkd^2d^i  (11) 
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with  (2)  and  (3).  Because  of  the  exchange  of  energy  between  the  particles  and  the  carrying  fluid,  this  ex¬ 
pression  of  does  not  provide  the  conservation  of  this  quantity.  If  we  consider,  in  addition  to  the  initial 
assumptions,  that  the  diameter  of  the  particles  is  very  small  compared  to  the  characteristic  length  (L)  of  the 
flow,  and  thanks  to  an  asymptotic  expansion  in  terms  of  the  small  parameter  D/L  «  1,  we  obtain  a  relation 
that  allows  to  prove  the  il-theorem. 


The  functions  'H{xi,t)  and  are  defined  by  the  relations  (12)  and  (13)  respectively  : 


=  J  f{^i,xi,t)ln  [/(^i,a:l,f)j  d^i 

(12) 

$H{x\,t)  =  J  ^if{^i,xi,t)ln  [/(6,ari,f) 

(13) 

Classically  the  relation  (14)  is  obtained. 

—n{xi,t)  -f  —  ■$H{xi,t)  =C{ln{f{ii,xi,t)) 

(14) 

If  we  consider,  in  addition  to  the  initial  assumptions,  that  the  diameter  of  the  particles  is  very  small  compared 
to  the  characteristic  length  (L)  of  the  flow,  and  thanks  to  an  asymptotic  development  of  (9)  in  terms  of  the 
small  parameter  D/L  «  1,  we  obtain  the  following  relation  : 

d  d  ->  —D^e^  ff  -  / 

-  /«■;,«)/(&, X-,))  In  I dU^dfi  +  0(0”)  (15) 

Easily,  we  proove  that  the  first  order  of  the  previous  expression  is  lower  than  zero.  It  means  that  the  quantity 
d'H/dt  +  (9./5aii).0/f  is  lower  than  zero  if  the  particles  are  sufficiently  small.  Consequently,  the  iJ-theorem  is 
true  for  small  particles  compared  to  the  characteristic  length  (L)  of  the  flow.  In  the  case  of  punctual  collisions, 
the  balance  law  of  the  %  quantity  is  reduced  to  the  first  order  of  the  previous  expression.  Then,  we  are  in  the 
same  conditions  as  [14]  but  we  assume  additionaly  that  the  collisions  are  weakly  inelastic.  It  allows  to  prove 
that  d'Hjdt  -I-  {d./dxi).(f)H  <  0  and  therefore  that  the  i?-theorem  is  true  in  that  case.  These  results  can  be 
generalized  to  the  case  of  a  bounded  medium  subject  to  the  choice  of  the  boundary  conditions. 


II  BGK  MODEL 

Because  of  the  complexity  of  the  non  linear  structure  of  the  collision  integral,  the  Boltzmann  equation  is 
very  difficult  to  solve  and  to  analyse.  It  is  therefore  interesting  to  introduce  a  simplified  model  equation  to 
study  qualitatively  its  solution.  There  are  principally  two  kinds  of  methods  to  obtain  these  models  from  the 
Boltzmann  equation.  The  first  consists  in  a  linearization  of  the  Boltzmann  equation.  We  find  this  approach  in 
[8]  for  instance.  The  second  is  based  on  the  unlinearized  Boltzmann  equation.  The  obtained  model  equation 
has  several  properties  of  the  original  operator,  but  it  is  not  strictly  derived  from  it.  We  find  this  approach  in 
[1],  [12],  [11]  for  instance.  An  approach  by  a  method  of  the  first  category  seems  to  be  difficult  to  apply  in  the 
case  of  not  punctual  collisions  of  inelastic  spheres,  studied  here.  Indeed,  the  use  of  a  method  of  linearization, 
requires,  for  the  construction  of  the  model,  the  conservation  of  mass,  momentum  and  energy  which  is  not 
guaranteed  by  the  original  model.  In  this  work,  an  approach  by  the  second  method,  which  has  not  the  same 
requirements  as  for  these  properties,  will  be  retained.  The  method  introduced  in  [12]  seems  to  be  difficult  to 
implement  in  the  case  studied  here.  Indeed,  it  requires  the  choice  of  a  model  distribution  function  /m  which 
allows  to  impose  constraints  on  model  and  which  is  the  limit  distribution  function  when  the  Knudsen  number 
Kn  of  the  flow  aims  towards  0.  The  choice  of  this  function  is  complicated  enough.  For  instance,  the  distribution 
function  of  ”13  moments”  does  not  agree,  because  it  does  not  allow  directly  to  express  the  second  member 
of  the  Boltzmann  equation  according  to  macroscopic  quantities.  Furthermore,  we  look  here  for  an  equation 
simplified  to  study  Knudsen’s  layer  near  a  wall  {Kn  <<  1).  The  use  of  this  method  would  mean  so  directly 
postulating  solution.  We  will  proceed  in  a  different  way.  Here,  the  second  method  is  used. 
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A  Splitting  of  the  Boltzmann  equation 


We  first  separate  the  collision  operator  of  the  Boltzmann  equation  in  two  parts.  This  step  is  very  classical 
(see  for  instance  [10]).  Using  the  modelling  (9)  and  introducing  the  functionals  Ji  and  J2,  the  Boltzmann 
equation  (8)  may  be  written  : 

If  +  ^  ^  ^  i  f)}  (16) 

with  : 

A(6,^i,i)  =  ^  [[  i9-k)lgoixi  +  ^k,t)f(^',xi,t)f{C'2,xi+Dk,t)\dkd^2  (17) 

e  JJ{g.k)>o  K  ^  ) 

and  : 

^  (g-k)  ^goixi  -  ^k,t)f{^2,xi  -  Dk,t)'^  dkd^2  (18) 


The  first  part  (called  Ji)  is  a  distribution  function  after  collision  while  the  second  {J2f{^i,Xi,t))  is  a  distri¬ 
bution  function  before  collision.  An  expression  of  these  two  functionnals  is  obtained  thanks  to  an  asymptotic 
expansion  in  terms  of  the  small  parameter  D/L  <<  1.  We  obtain  : 


Ji{ii,xi,t)  =  D^go{xi,t)  [[  ^  {g-k)f' f'2dkd(2  + D^go{xi,t)  [f  ^  {g  ■  k)k  ■  /  ^dkd^2 

JJ(g-k)>0  JJ  {g-k)>0  OX\ 

+  .  ff  (^ .  k)f'f2kdkd{2  +  0{D^)  (19) 

OXl  JJ(3.k)>0 


and  : 


J2iii,xi,t)  =  D‘^goixi,t)  [[  {g-k)f2dkd^2- D^goixi,t)  [[  {g  ■  k)k  ■  ^dkd{2 

JJ(gk)>Q  JJ(3k)>o  0x1 

-  .  [[  (^ .  k)f2kdkd{2  +  0{D^)  (20) 

0X1  JJ(Qk)>0 


Firstly,  we  work  on  the  expression  of  J2{^i,xi,t)  . 


B  Modelling  of  ^2(^15^!,^) 

At  this  stage  of  the  calculation,  within  the  framework  of  the  not  linear  models  of  the  Boltzmann  equation, 
assumptions  are  usually  made  to  eliminate  g  of  the  integrals.  It  is  case  for  instance,  for  the  model  BGK  (  [1], 
[10]).  Indeed,  this  model  is  based  on  the  use  of  a  maxwellian  potential  which  allows  to  modify  slightly  the 
writing  of  J2.  Introducing  then  the  concept  of  “pseudo-maxwellian”  particles,  the  authors  make  indirectly 
hypotheses  on  g.  Then,  this  quantity  disappears  of  the  integrals  and  is  absorbed  by  means  of  a  constant 
frequency  collision.  In  the  case  studied  here,  this  kind  of  hypothesis  is  impossible  to  use.  Indeed,  the  dynamics 
of  the  inelastic  collision  of  two  particles  is  not  comparable  with  that  of  two  molecules  where  attraction  / 
repulsion  forces  play  a  dominating  role.  The  same  objections  can  be  set  against  the  use  of  a  Lennard-Jones 
potential,  often  used  in  the  kinetic  theory  of  gas.  In  order  to  be  consistent  with  the  assumptions  related  to  the 
collisions  and  supposing  that  particles  are  big  enough  to  be  able  to  neglect  the  interaction  potential  between 
particles  and  molecules  of  the  carying  fluid  presented  in  [13],  we  can  not  use  the  pseudo-maxwellian  particles 
assumption.  Besides,  with  a  rigid  spheres  potential,  which  is  the  most  consistent  with  the  retained  collision 
modelling,  the  frequency  of  collision  varies  with  the  particles  velocities,  and  this  variation  is  expected  to  be 
important  when  the  particles  velocities  are  high  (  [3]).  Objective  is  so  to  introduce  in  the  model,  a  collision 
frequency  depending  on  the  particles  velocities.  By  realizing  integration  in  fc  in  J2,  we  obtain  : 
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|fl(o(fi,i)7r  j  gf2d^2  -  ^^9o{xi,t)  j  g- 


df2  2ttD  dgo{xi,t) 


6  5x1 


■  I  gf2d{2^+OiD^) 
(21) 


with  <?  =  1 1  5  1 1  ■  In  the  aim  to  obtain  a  simplified  expression  for  J2 ,  we  use  a  similar  method  to  that  of  the  usual 
case.  We  assume  that  J2(^i,xi,f)  may  be  written  under  the  form  ; 


J2{ii,xi,t)  =  i^goixi,t)ire{ii)^i  j  f2d^2 

2^Ddgo{x^,t)  I 


^^go{Xl,t)e{^l)il  ■  ^  /  f2d^2 


f2d^2\+0{D*)  with  6=IUi 


and  then 


6  9xi 


J2i^i,Xi,t)  =  \  goixi,t)Tre{ii)^in{xi,t)  - 


2ttD  5go(xi,f) 

6  dxi 


,*)} 


•0(6)ein(xi,f)  +0(i?4) 


That  is,  always  with  D  «  L  : 


J2  =  TT  1 9o{xi  -  -  ^)  +  0{D^) 


(22) 


(23) 


(24) 


Then  we  set  iV(xi  —  —^)  =  go{xi  —  ^-)n{xi  —  and  taking  this  into  account,  J2  is  modelled  by  : 

?i  241  4i 

J2  =  i'{xi-^p-,t)  with:  v{xi-^^,t)  =  -K  D'^  go{xi-^-,t)n{xi-‘^^^p-,t)  (25) 

4i  4i  d4i  341 

where  n  is  the  number  density  and  0  a  function  of  the  particle  velocity. 


C  Full  operator  modelling 

In  accordance  with  the  previous  section,  the  new  equation  is  : 

As  in  the  derivation  of  the  usual  BGK  model,  we  put  : 

Jl  =  \  l^{Xl  +  fo{^l,Xi,t) 

4i 


(26) 


(27) 


where  /o(4i)®i,t)  is  a  Maxwellian  distribution  function.  This  choice  expresses  the  tendency  of  the  suspension 
to  return  towards  an  equilibrium  state.  It  is  confirmed  by  the  existence  of  an  H  theorem  for  the  original 
Boltzmann  equation.  To  satisfy  in  best  the  properties  of  the  original  collision  operator,  it  is  advisable  to  bring 
the  biggest  care  to  the  choice  of  this  function.  It  is  chosen  as  the  maxwellian  satisfying  the  following  three 
integral  properties: 


l>  Two  properties  allowing  the  conservation  of  mass  and  momentum  : 


Xi  + 


D^i 


'ii 

with  tp  =  ^1  and  ip  =  1. 


t)fo{^i,xi,t)-^d^i  =  f  e^z/(4i,xi  -  I^,t)/(4i,xi,t)V'd4i 
Jfi  4i 


(28) 
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[>  One  assumption  about  the  energy  balance  : 


/ 

Jii 


[ 


(29) 


The  density,  velocity  and  energy  which  appear  in  are  not  the  five  local  moments,  but  some  fictitious 

local  parameters  related  to  the  five  moments  of  /  weighted  with  1/(^1, The  obtained  BGK  model  is 
therefore  : 


df  df  ^  df 


D^i 


<  i'{^i,xi  +  —,t)fo{^i,xi,t)-e  u{^i,xi  -  —,t)f{^i,xi,t) 


(30) 


This  model  is  similar  to  the  usual  BGK  model  in  the  case  of  non  punctual  and  inelastic  collisions  of  identical 
hard  spheres. 


One  can  show  that  by  adding  some  assumptions  to  this  model,  one  finds  known  models.  Supposing  that 
collisisons  are  punctual  and  elastic  (e  =  1)  and  supposing  that  the  shield  effect  is  negligible  {go  =  1),  the 
equation  (30)  may  be  written  : 

+  =  W6,^i,li)  [fo{ii,xi,t)  -  /(^i,xi,t)|  (31) 


with  : 


^  yTT  /  j  ^ 

t^{^i,xi,t)  =  —^^i6{^i)n{xi,t)  (32) 

It  is  a  generalization  of  the  BGK  model  with  a  collision  frequency  which  depends  on  the  particles  velocity  ( 
[3],  [11]).  The  previous  relation  is  a  particular  expression  of  the  collision  frequency  v. 

Furthermore,  if  the  assumption  of  pseudo-maxwellian  particles  is  made,  the  collision  frequency  does  not 
depend  any  more  on  the  particles  velocity  and  it  mays  be  written  : 

z/(|*i  ,xi,t)  =  A  n{xi ,  t)  (33) 

where  the  value  of  A  may  be  found  in  [10]  for  instance.  The  obtained  model  is  then  the  usual  BGK  model. 


Ill  PROPERTIES  OF  THE  BGK  MODEL  EQUATION 


A  Equation  of  moments 


In  a  very  classical  way  in  kinetic  theory,  the  equation  of  moments  is  written  : 

ipd^i  =  C{tp) 


df  z-  df 
dx 


1  9a/ 


with  : 


‘"W  =  ^  j  +  ^-,t)fo{ii,Xi,t) -e^i'{ii,xi  -  ^,t)/(a,fi,t)|  V’dfi 


(34) 


(35) 


Due  to  the  assumptions  (28),  the  conservation  of  mass  (V’  =  1)  and  of  momentum  {tp  =  ^1)  are  automatically 
verified.  On  the  other  hand,  as  for  the  original  operator,  this  equation  does  not  give  the  conservation  of  energy 
{tp  =  a).  Model  supplies  so,  at  this  level,  the  same  type  of  properties  as  the  original  equation  except  for  the 
momentum.  However,  if  the  assumption  (which  we  interhook  mainly  here)  of  small  particles  is  made,  the  BGK 
equation  and  the  Boltzmann  equation  have  the  same  properties.  Under  this  assumption,  it  remains: 
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(36) 


CW 


+  ^  J  ^  (^fo{ii,xi,t)  +  e‘^f{^i,xi,t)^  dil  +  OiD"^) 


B  H  Theorem 

The  functions  H{xi,t)  and  (j)H{x\,t)  are  defined  in  the  same  way  as  in  the  first  section  by  the  relations  (12) 
and  (13)  respectively.  In  the  same  way  as  in  the  first  section,  the  following  law  of  balance  for  the  %  quantity 
is  obtained  : 


d 


d 


1 


—n{xi,t)  +  -^-(i)H{xi,t)  =  ^  /  <  1/(^1, fi  + -—,t)/o(^i,fi,t) 


dt 


dxi 


a 


e  ^  [l  +  /n(/)]d^i 


(37) 


In  the  case  of  small  particles  and  using  the  conservation  of  mass  ,  it  remains  : 

^'H{xi,t)  +  ^  •  $H{x1,t)  =  ^  j  [Hii,xi,t)foi^i,xi,t) -e^iy{ii,xi,t)f{ii,xi,t)'^ln  d^i 

+  ^j  {^(.ii^xi,t)fo{ii,xi,t) -e^i'{ii,xi,t)f{ii,xi,t)'^ln{fo)d^i 

+  0(0)  (38) 

In  the  cas  of  weakly  inelastic  collisions,  we  set  e  =  1  —  ry  with  r}  «  1.  Consequently,  the  previous  expression 
is  reduced  to  ; 

■  ^Hix'\,t)  =  j  |t/(6,fi,t)/o(|'i,fi,t)  -  J^(6,ri,<)/(|'i,xi,t)|  In  d^i 

+  J  [i^{(i,xi,t)fo{ii,xi,t)-u{(i,xi,t)f{^i,xi,t)^ln{fo)d^i 

+  0{D)  +  0{ri)  (39) 

The  function  fo  being  maxwellian,  ln{fo)  may  be  written  according  to  the  quantites  1,  and  Consequently, 
the  second  inegral  of  the  previous  expression  may  be  written: 

j  {j^(6,xi,t)/o(6,xi,t)  -  n{^i,xi,t)f{ii,xi,t)^ln{fo)d^i 

=  J  {i"(^i,Xi,t)/o(^i,Xi,t) -i/(|*i,Xi,t)/(fi,:Ei,t)|  [ao  +  dl -fi +0:26^]  (40) 

By  using  the  assumption  (29),  it  remains,  in  the  second  integral  of  the  previous  expression,  only  the  terms  ”1” 
and  which  nullify  by  using  a  linearization  of  the  first  order  of  relations  (28): 

J  |i^(fi,xi,t)/o(^i,xi,t)  -  i/(fi,xi,t)/(fi,xi,t)|  in(/o)dfi  =0  +  O{D)  (41) 

Consequently,  the  balance  equation  for  7i{xi,t)  is  reduced  to  : 


<Aff(xl,t)  = 


J 


Xl,t)foi^l,Xi,t)  -  Z/(^i,Xi,t)/(^l,Xi 


(s) 


d^i 


+  0{D)  +  0{ri) 


(42) 
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The  collision  frequency  is  necessarily  greater  than  0  and  the  term  : 


is  always  lower  than  0.  We  have  therefore  at  the  first  order  : 

d  d 

—nixi,t)  +  —  ■  <  0  (44) 

Consequently,  if  we  consider  small  particles  and  weakly  inelastic  collisions,  we  have  an  H  theorem  for  the  BGK 
model,  that  is  the  same  result  as  for  the  original  Boltzmann  equation. 


IV  CONCLUSION 

The  problem  of  a  two-phase  dispersed  medium  is  studied  within  the  kinetic  theory.  In  the  case  of  small 
and  undeformable  spheres  that  have  all  the  same  radius,  a  model  for  the  collision  operator  of  the  Boltzmann 
equation  is  proposed.  The  collisions  axe  supposed  instantaneous,  binary  and  inelastic.  The  obtained  collision 
operator  allows,  to  prove  the  existence  of  an  H  theorem  in  the  case  of  a  suspension  of  small  hard  spheres  having 
binaray,  instantaneous  and  weakly  inelastic  collisions.  Because  of  the  complexity  of  the  non  linear  structure 
of  the  collision  integral  of  the  Boltzmann  equation,  it  is  interesting  to  introduce  a  BGK  model  equation.  All 
the  basic  formal  properties  of  the  original  operator,  are  retained  by  this  model  equation  in  the  case  of  small 
particles.  Consequently,  this  model  is  consistent  with  the  original  kinetic  description  and  it  will  be  very  usefull 
in  the  study  of  some  complicated  problems  like  the  Knudsen  layer  for  instance. 
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